A quaternionic generalisation of the Riccati 
differential equation 



Viktor Kravchenko, 

Univ. do Algarve, Dep. Math., Campus de Gambelas, 8000 Faro, Portugal 

Vladislav Kravchenko, 

Esc. Sup. de Ing. Mec. y Elec. del Inst. Polit. Nac., Dept. de Telecom., 
Unidad Zacatenco, CP. 07738, D.F., Mexico 

Benjamin Williams 

Esc. Sup. de Fis. Mat. del Inst. Polit. Nac, 
Unidad Zacatenco, CP. 07738, D.F., Mexico 

February 5, 2008 

Abstract 

A quaternionic partial differential equation is shown to be a gener- 
alisation of the traditional Riccati equation and its relationship with the 
Schrodinger equation is established. Various approaches to the problem of 
finding particular solutions to this equation are explored, and the generali- 
sations of two theorems of Euler on the Riccati equation, which correspond 
to this partial differential equation, are stated and proved. 

1 Introduction 

The Riccati equation 

du = pu + qu + r, 

where p, q and r are functions, has received a great deal of attention since a 
particular version was first studied by Count Riccati in 1724, owing to both its 
peculiar properties and the wide range of applications in which it appears. For 
a survey of the history and classical results on this equation, see for example 
ID, § and Q. This equation can be reduced to its canonical form W, 

dy + y 2 = -v, (1) 



and this is the form that we will consider. 



One of the reasons for which the Riccati equation has so many applications is 
that it is related to the general second order homogeneous differential equation. 
In particular, the one-dimensional Schrodinger equation 

-d 2 u — vu = (2) 

where v is a function, is related to the (|l|) by the easily inverted substitution 
y = £p This substitution, which as its most spectacular application reduces 
Burger's equation to the standard one-dimensional heat equation, is the basis 
of the well-developed theory of logarithmic derivatives for the integration of 
nonlinear differential equations [ pd| . A generalisation of this substitution will 
be used in this work. 

A second relation between the one-dimensional Schrodinger equation and 
the Riccati equation is as follows. The one-dimensional Schrodinger operator 
can be factorised in the form 

~d 2 ~v(x) = -(d + y(x))(d-y(x)) 

if and only if (]]]) holds. 

Among the peculiar properties of the Riccati equation stand out two theo- 
rems of Euler, dating from 1760. The first of these jl3) states that if a particular 
solution yo of the Riccati equatuion is known, the substitution y = yo + z re- 
duces (|l|) to a Bernoulli equation which in turn is reduced by the substitution 
z — — to a first order linear equation. Thus given a particular solution of the 
Riccati equation, the general solution can be found in two integrations. The 
second of these theorems Q states that given two particular solutions yo, y\ 
of the Riccati equation, the general solution can be found in the form 

y _ kyoexpjfyo - y x ) - y x 
kexp(J yo-yi) - 1 

where k is a constant. That is, given two particular solutions of (Q), the general 
solution can be found in one integration. 

Other interesting properties are those discovered by Picard and Weyr Jl3| | . 
The first is that given a third particular solution 7/3, the general solution can be 
found without integrating. That is, an explicit combination of three particular 
solutions gives the general solution. The second is that given a fourth particular 
solutions 2/4, the cross ratio 

(2/1 -2/2X2/3-2/4) 

(2/i - 2/4)0/3 - 2/2) 

is a constant. 

This article is concerned with a quaternionic generalisation of the Riccati 
equation and versions of the above-mentioned theorems of Euler corresponding 
to this generalisation. Some necessary notation will be introduced in Section ??. 
In Section [3] we propose the quaternionic generalisation of the Riccati equation, 
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which is shown to be a good generalisation for various reasons, including that 
it is related to the three-dimensional Schrodinger equation 



Au + vu = (4) 

(where A is the three-dimensional laplacian) in the same ways as the Riccati 
equation is related to (0) . In Section |], we turn our attention to cases in which 
particular solutions of (g) can be found, some of which differ considerably from 
the one-dimensional case. Finally, in Section ||, generalisations of Euler's 
theorems will be stated and proved. 



2 Preliminaries 

The complex numbers and complex quaterions are denoted by C, H(C) respec- 
tively. The latter consists of elements of the form 

3 

a = ^2 a k ik 

k=0 

where the 8* eC and the base units ik satisfy the following rules of multiplica- 
tion 



i Q = i = -ik, ioik = ioik = ik, k= 1,2, 3, 
iii 2 = -i 2 h = «3, «2«3 = -i3«2 = h, «3«i = -hh = 

The complex unit i commutes with the ik- Frequently it is useful to consider a 
quaternion a as being the sum of a scalar and a vector part, denoted respectively 

3 

ao := Sc(a), it := Vec(a) = cikik- 

k=l 

Conjugation is defined as follows, 

a := ao — a 

and the modulus is 

| a. | 2 = a ■ a = dp + af + a\ + a\ 

so that in particular it 2 = — |lz*| 2 . 

Note that in terms of scalars and vectors, the quaternionic product can be 
written 

ab = (ao + ~c?)(bo + b ) = ao&o + aa b + bolt — (it, b \ + it x 6 
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where (a, b) is the standard inner product and [a x b] the standard cross product 
in R 3 . In particular 

{o>,T} = -2(a>,T) (5) 

where {a, b} = ab + ba is the standard anticommutator. 

In what follows functions g : £1 — » H(C) will be considered, where f2 is some 
domain in R 3 . The Moisil-Theodoresco operator D is defined on differentiable 
functions g as follows: 

3 

Dg = ^] ikdkg 

k=l 

where dk — g^-. Due to properties of the quaternionic product, this can be 
written 

Dg = — div if + grad go + r °t ~ff ■ 

Thus it follows that 

Sc(zrr) = -div 

Vec{D~(f) = rotlf, 

and for scalar functions u 

Du = gradit. 

The theorem of Leibnitz for this operator is the following: given a differen- 
tiable scalar function u and a differentiable quaternionic function g, 

D(ug) = D(u)g + uD{g). (6) 

The logarithmic derivative (in the sense of Marchenko) of a scalar function 
u such that u ^ in Q, is defined as 

du = u^Du. 

The function du is a vector. The derivative is logarithmic in the following sense: 
given two scalar functions ui,i«2 that do not vanish in f2, formula (|^) implies 
that 

d{u\U2) = du\ + Bu2- 
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3 A three-dimensional generalisation of the Ric- 



cati equation 

The substitution of the Jackiw-Nohl-Rebbi-'t Hooft ansatz [0 in the self-duality 
equation can be written H in the following quaternionic form 



d t g + Dg + \g\ 2 = (7) 

where the subscript t denotes differentiation with respect to time. This equation 
has obvious formal similarities with equation (0). In || the relation of (0) to 
the Filter operator 

d t + D 

was shown. In particular, it was shown that for any / € ker(dt+D), the function 

2(grad/ - div /*) 
fo 

is a solution of (0), that is, a class of instantons was obtained. In what follows we 
will concentrate on the case of time independent, purely vectorial quaternionic 
functions, but the nonhomogencous equation will be considered. 

The following result generalises to three dimensions the relation, mentioned 
in the introduction, between the one-dimensional Schrodinger operator and the 
Riccati differential equation via the logarithmic derivative. 

Proposition 1 ip is a solution of fy) if and only if f := dip is a solution of 

Dj + J 2 = v (8) 

Proof. Suppose that there exists a function tp such that 7 = dip. Applying 

(§) gives 

D7 =\ (Vip,S7ip) ~ -Aip, 

so that —^Aip — v, or equivalently Aip +vip = 0. Conversely, given a solution 
tp of (f|) , / = dip is a solution of (§) . ■ 

In H and S it was shown that the three-dimensional Schrodinger operator 
can be factorised in the following way 



-A-vI= (D + M*){D -M f ) 
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where / is the identity operator and M f q := q f , if and only if equation (0) 
holds. 

Thus the two relations between the Riccati equation ([l]) and the one-dimensional 
Schrodinger equation (Q) mentioned in the introduction have natural counter- 
parts relating (||) and the three-dimensional Schrodinger equation (Q). These 
relationships suggest that (||) can be considered a good generalisation of (jl]) . It 
should also be noted that if / = fk{xk)iki then (|J) is reduced to 

dkh + fk = 

That is, a one-dimensional solution of (|8j) is a solution of (|l|). Equation (Q) will 
be refered to as the Riccati PDE. 

Note that the scalar and vector components of equation (||) are respectively 

-divT + 7 2 = v, 

rot 7 = 0. 

The second equation implies that for a simply-connected domain f2, there exists 
a scalar function ip such that / = grad ip. Substituting this in the first equation 
gives 

Atp+ (V<^, Vp) = -v. (9) 

This equivalence of the Riccati PDE with a scalar elliptic partial differential 
equation will be used frequently in what follows. 

It should also be noted that if the function v in (||) is zero, the substitution 
/ = dtp reduces the equation to 

= 0. 

Thus the homogeneous Riccati equation can be solved explicitly, and its solu- 
tions are of the form dip, where ip € ker A. This generalises the highly restricted 
class of solutions of the homogeneous equation (|IJ), which are precisely functions 
of the form 

1 

x + c' 

c constant. 



4 Particular solutions of the equation 



In the next section it will be shown that given one solution of the Riccati PDE 
it can be linearised. Thus in this section we discuss some possibilities for 
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obtaining particular solutions of (||) . First we note that if the function v is of 
the form 

v(x) = vi(xi) + v 2 (x 2 ) + v 3 (x 3 ), 

then assuming that / = fi(xi)ii + ,f 2 {x2)i2 + /3(^3)*3, equation (||) reduces to 
the system of Riccati ordinary differential equations 

dkfk + fk=-v k , k = 1,2,3. 

Thus in this particular solution of (|§|) can be found if and only if each 

of the above equations can be solved. Obviously, if any of the Vk's are zero, 
this task is greatly simplified. This situation corresponds to the Schrodinger 
equation ^ with potential v of the form given above, in which case the variables 
can be separated. 

The existence of a large class of solutions of the homogeneous equation, as 
described in Section [| motivates the following procedure, which reduces (||) 
to various scalar differential equations. Substituting the sum of two functions 
/ 1 ! / 2 € C 1 (fl) into equation (||) gives the following possible decomposition: 

d7 1 + 7 2 1 = v 1 , 

{ f l, ^2} = v 3 , 

v = vx + v 2 + Vs. 

In particular, if / 1, / 2 are solutions of the homogeneous equation, this is 
reduced to 

{Ti,7 2 } = -2(dtpi,8ip 2 ) = v, (10) 
where ip±, ip 2 are harmonic functions. 
If tpi = xi, / 1 = ^l, this becomes 

di<p 2 = -\^x X (p 2 



which has solution 



ip2 = A(x 2 , £3) exp(— - J vxidxi 



where A(x 2 ,x 3 ) is an arbitrary function. Thus if if2 G ker A and if2 satisfies 
the above equation, the sum 

— + dif 2 

Xl 

is a particular solution of the Riccati PDE. 
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If instead of choosing ipi as above ip\ — ip2 is substituted in (TO), the eikonal 
equation 

2(cV) 2 = -v (11) 

results. Thus for a scalar function ip e ker A3, which is also a solution of the 
above eikonal equation, / = 2dtp is a solution of (||) ■ 

Example 2 Let tp be the fundamental solution of the laplacian A, 

1 



47T \ x\ 



This function is harmonic and positive in any domain Q which does not include 
the origin. Furthermore it satisfies equation Mm with 



1 

v = — 



2 • 

X 



Thus 



Ixl 2 



is a solution of , 



5 Generalisations of Euler's theorems on the Ric- 
cati equation 

We now state and prove the generalisations to the Riccati PDE of Euler's the- 
orems on the Riccati equation that were mentioned in the introduction. The 
first of these theorems states that given a particular solution of the Riccati 
differential equation, the equation can be linearised. 

Proposition 3 (Generalisation of Euler's first theorem) Let h — grad£ be an 

arbitrary particular solution of (Qj. Then 

7 = l? + ~ti (12) 
is also a solution of where if = and ^ is a solution of the equation 

A* + 2 (V£, W) = 0, (13) 

or equivalently of 

div(e 2< = W) = 0. (14) 
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Proof. Substituting (|T|) in (g) gives 

D If + {hi ~g>} + Y 2 = 

or alternatively, using (g) 

D~g>-2(t, 1?) + T 2 = 0. (15) 
Note that, as in (g), the vector part of ( [l5| ) is rot if = 0, so that 

~c? = grad<& 

for some function $. If W = e*, this is equivalent to 

T = <9*- 

Equation ([l5|), written in terms of \&, is 



so that (15) is equivalent to 

A* + 2 (V£, W) = 0. 

Noting that 

div(e 2 «W) = 2e 2 « (V£, V*) + e 2 « A* = e 2 «(A* + 2 (Vf , V*)), 
this equation can be rewritten in the form 

div(e 2S W) = 0. 



Equation ( |l3| ) is the well-known transport equation, which appears for ex- 
ample in the Ray Method of approximations of solutions to the wave equation, 
coupled with the eikonal t 
tions, for example in electrostatics 
and is the electric field potential, and as the continuity equation of hydrome- 
chanics in the case of a steady flow, where e 2 ^ is the density of the medium 



equation jl|]. Equation (O) appears in various applica- 
ctrostatics |lQ |, where e 2 Ms the dielectric permeability 



Remark 4 From (l^ ) we have that 

e^V* = rot V 
for some vector-valued function ~s*, or 

V<3> = e" 2? rot ~s*, 
where ~s* must satisfy the condition 

rot(e~ 25 rot^*) = 
as e~~ 2 ^ rot if must be the gradient of some function. 
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As mentioned in the introduction, given two particular solutions y\ 1 1/2 of 



the Riccati ordinary differential equation the general solution can be found in 
one integration. A natural question is whether this property extends to the 
Riccati PDE. The form (3) of the general solution found in this case, suggests 
a similar substitution in (3). This line of reasoning gives the following result. 

Proposition 5 (Generalisation of Euler's second theorem) Let h 1 = grad£i 

and h 2 = grad^2 be two particular solutions of ^). Then there exists a scalar 
function if such that 

7j\w - 7^2 

/ = V(p = 



w — 1 

is also a solution of where w = Ae^ 1 ^^ 2 , A G 
Proof. Equation <M) is equivalent to 



Atp+ (Vtp, Vtp) = -v (16) 

where ~f = gradt^. The scalar functions £1 and £2 are two solutions of (|l6|). 
Substituting the expression 

V ^ - V ^ 2 (17) 
w — 1 

for V</? into this equation, where w is a scalar function, gives 

a v = div (Y^lk 

\ w — 1 

_ (w - l)(wA£i + Vtu ■ V^i - A£ 2 ) - Viu • (wV£i - V£ 2 ) 
~ (w - l) 2 

and 



(iy-1) 2 



Simplifying and using the fact that £1 and £2 are solutions of ([16]), the equation 
is reduced to 

Vlogw = — = V(£i -6), 

so that 
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for an arbitrary constant A £ C. It remains to show that the expression ( p"7| ) 
is the gradient of some scalar function <p. This is the case if the rotational 
of (p"7| ) disappears. This is shown as follows, where the identities rotVt^ = 
0, / x / = are used. 

rot ,v^-v6U ro Vv^A ro /v6 



w — 1 / \ w — 1 / \w- 1 

grad x V£i H rotV£i 

...i w — 1 



w — 


1 


w 




w — 


1 


1 




w — 


1 


w 





grad — ^— x V£ 2 — r rot V£ 2 

w — 1 

grad — — - x VCi - grad — — x V£ 2 
w — 1 w — 1 

/ tt^ x + - tt^ x V6 

(w — i) z (w — \y 

(u> — \y 

~{W- l) 2 

0. 



V(6^i)xV(6-a) 



It must be noted that the new solution gained is not necessarily the general 
solution, but a larger class of solutions. The following example illustrates this. 

Example 6 Consider ^) with v = — 1, 

d7 + 7 2 = -i. 

Two solutions of this equation are h\ = i\ = gradxi, /i 2 = z 2 = gradaja- 
Applying the above result gives the class of solutions 

u Ae Xl ~ X2 — in 

-^P AeC, 

Ae Xl - X2 - 1 

however a third solution h% — 13 is not included in the above expression as a 
special case. 
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